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Analyzing simulated EAS events generated with the CORSIKA code, the characteristics of lat- 
eral distribution of electrons in EAS around the knee energy region of the primary energy spectrum 
have been studied and compared with experimental observations. The differences between the EGS4 
and the NKG output of CORSIKA in respect to electron radial density distribution have been in- 
vestigated. The relation between lateral and longitudinal age parameters has been studied after 
introducing the notion of the local age parameter that reflects the profile of the lateral distribu- 
tion of electrons in EAS. The present analysis motivates the inclusion of the lateral shower age 
in a multiparameter study of EAS to extract information on hadronic interactions and primary 
composition. 
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I. THE DIFFERENT AGE PARAMETERS 

The concept of shower age was introduced in cas- 
cade theory to describe the stage of development of 
an electromagnetic (e.m.) cascade. A synthesis sum- 
marizing the works of Greisen and Nishimura-Kamata 
under Approximation B of cascade development [l[ in 
respect to shower age is the following: the longitudinal 
age sl is defined here as 



sl 



3t 



t + 2ln(E/eo) ' 



(1) 



where Eq is the energy of the primary photon gener- 
ating the cascade, t is the atmospheric (divided by the 
electron radiation lenth in air taken as 37.1 g-cm~ 2 ), 
eo being the critical energy of 82 MeV. 

In this theoretical context the lateral density distri- 
bution of cascade particles given by Nishimura and 
Kamata can be approximated by the well known 
Nishimura-Kamata-Greisen (NKG) structure func- 
tion, 



f(r) = C(s ± )(r/r m y±- 2 (l+r/r m ) s ± 



-4.5 



(2) 



where r is the radial distance measured from the 
EAS core, r m is the Moliere radius, s± is the lateral 
age : the normalization factor C(s_i_) is given by 



C(s ± ) 



r(4.5 - 8 ±) 



27rr(s_L)r(4.5 - 2s j 



(3) 



implying that for the density pNKc{ r ) = N e f(r) 
thanks to the properties of the Eulerian function. 

The relation sl = s± was initially considered to 
hold for pure e.m. showers and it was admitted that 



the average steepness of the profile of the lateral dis- 
tribution has a correlation with the longitudinal de- 
velopment. 

Later the 3D diffusion equations were solved again 
by Uchaikin and Lagutin using adjoint equations and 
an improvement of NKG function was proposed by 
modulating r m to sl as follows: 



Pei(r) = (mr m ) 2 p N KG(r/m) 



(4) 



with m = 0.78 - 0.21s L . 

The validity of this approach was demonstrated Q 
by pointing out that the experimental distributions 
in EAS are steeper than the Pnkg but are in better 
agreement with Monte Carlo calculations of Hillas. 

For better estimation of shower age from the ex- 
perimental distributions, one of us (Capdevielle) in- 
troduced the notion of the local age parameter (LAP) 
Shoe 'M- From two neighbouring points, i and j, one 
can give a lateral age parameter for any distribution 
f{x) (where x = that characterises the best fit by 
a NKG- type function in [xi,Xj] : 



HF t] Xf 3 Y^) 
kiiXijYij) 



(5) 



where i^y = f(n)/f(rj), Xij^n/rj, and 
Y i j = (x i +l)/(xj+l). More generally, if — > r 
this suggests the definition of the LAP s(x) (or s(r)) 
at each point : 



3 ' 



s(r) = 



1 



2x + 1 



, <91n f 

(x + l)g^ + (2 + n )x + 2) (6) 



If /3o=4.5, fNKa{ r ) w ith s=s(r) can be used to fit 
/ in the neighbourhood of r. 
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The identification s(r) — Sij for r = r ' 2 rj remains 
valid for the experimental distributions as long as they 
are approximated by monotonic decreasing functions 
versus distance: A typical behaviour of s(r) was in- 
fered with a characterised minimum value of the pa- 
rameter near 30 — 50 m from the axis followed by a 
general increase at large distance and it suggests a 
relation sl ~ 1.25 — 1.3sj_. 

After verification of the behaviour of the LAP 
through experimentally observed lateral distributions 
and particularly a detail study of the parameter with 
the Akeno data 0] , this approach was validated by the 
rapporteurs of the ICRC from 1981 to 1985 [jj. This 
procedure was also used in extension to calculate the 
radio effect of very large EAS @ . 
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II. 3D SIMULATION OF EAS 

In the present work, the high energy (above 
80GeV/n) hadronic interaction model QGSJET 01 
version lc has been used in combination with the low 
energy (below 80GeV/n) hadronic interaction model 
GHEISHA (version 2002d) in the framework of the 
CORSIKA Monte Carlo program version 6.600 Q to 
generate EAS events. The simulated shower library 
produced consists mainly of 10, 000 EAS for each pri- 
mary species Proton, Helium and Iron in the primary 
energy interval of 10 14 eV to 3 x 10 16 eV. 

Taking the opportunity of calculating simultane- 
ously the e.m. component via both the options, the 
EGS and the NKG option (implemented following [3[ 
by exploiting the relation (4) in the subcascades treat- 
ment) as admitted in CORSIKA package, we have car- 
ried out parallel simulations. We ascertain that when 
the calculation of the electron component is carried 
out with relation (4), the situation with the NKG in- 
spired procedure implemented in CORSIKA is more 
close to the experimental data and also to the cal- 
culation with the EGS, as shown in fig.l. The NKG 
option gives a slighly larger density with steeper radial 
distribution in compare to the EGS option. A small 
density excess appears for pure electromagnetic cas- 
cades near the axis for the Corsika NKG option. Such 
an excess appears also in proton induced air show- 
ers. However, over a large band of densities between 
the radial distance 10m up to 100m from the axis, a 
tolerable agreement is noticed. 

For Fe primaries, both the options produce older 
profile near the axis with density excess (particularly 
in the case of NKG option) between 2 — 10m distance. 
The average energy of the positrons is quite small in 
the case of iron initiated showers and the cross section 
of positron annihilation becomes more important for 
the lower part of the cascade. This effect is probably 
enhanced by the longer path of the electrons in the ge- 
omagnetic field (and larger energy loss by ionisation) . 
Therefore the NKG option is not a very good option 



FIG. 1: Comparison of radial distribution of electrons for 
CORSIKA NKG and EGS options for proton and Fe pri- 
maries. 



for the simulation of heavy nuclei initiated inclined 
showers after the maximum development of cascade. 
Near vertical showers in low altitudes are nevertheless 
acceptable. Conversely for very high energy primaries 
the average energy of the positrons remains impor- 
tant and for proton and photon showers NKG option 
is still useful to calculate a large number of cascades 
in a short time. 

At larger distances a slight deficit in densities ap- 
pears with the NKG option; it seems to come from the 
different treatment of the multiple coulomb scattering 
in the EGS and also may due to the enhanced path 
of the muons in the geomagnetic field which results in 
greater loss of energy by ionisation and subsequently 
their decay produces more electrons. The EGS takes 
into account the photoproduction inside the e.m. sub- 
cascades. 

Through a smaller sample of events simulated in 
the energy band 10 18 — 10 20 eV [§(, we have also ob- 
served that the lateral distributions calculated at dis- 
tances lower than 300m from the axis with the NKG 
option as well as the total longitudinal development 
do not differ much from those obtained with the EGS. 
Consequently, the NKG option in CORSIKA remains 
useful for faster calculations at ultra high energy initi- 
ated by a nucleon or nuclei. This circumstance allows 
the calculation of radio effect as in [6j as well as the 
fluorescent component via the NKG option. The Lan- 
dau Pomeranchuk Migdal effect (LPM)(not included 
in the NKG option) limits however the employment 
of the option beyond 10 18 eV for 7 primaries. 

The fluctuations in lateral shower age are much 
larger for proton initiated showers compared to those 
initiated by heavier primary. 

This distribution (fig. 2) can be fitted by an Extreme 



Cll 



XVI International Symposium on Very High Energy Cosmic Ray Interactions 
ISVHECRI 2010, Batavia, II, USA (28 June 2 July 2010) 



3 



s fluctuations and EVD distribution Electron latera distributor 




s 



FIG. 2: fluctuation of shower age at fixed shower size 7V e = 
10 4 fitted by an E.V.D. function 4>{s) (vertical showers, sea 
level, p primaries) 

Value Distribution (E.V.D) denned through 

cj>( S ) = - exp(±^^ - e ±Ll ^ 1 ) (7) 
a a 

where the parameters [x and a are related to the 
average size s and its variance V s by s = fj, ± 0.577 a 
and V s = 1.645a 2 (in the case of the histogram of 
fig.2, s = 1.495 and a = 0.07). 

Comparing the results of the EGS and the 
NKG options with Kascade data @ in the bin 
Log(N e ) [3.9, 4.3], we observe that there are small dis- 
crepancies in densities even with the Lagutin NKG 
modified formula taking s = 1.404, N e = 11829 as 
shown on fig. 3 

Those mentioned parameters (for Lagutin formula 
(4)) corresponds respectively to the average lateral age 
parameter and the average size in the size bin consid- 
ered (calculated with the size spectrum in Karlsruhe). 

III. LAP AND N e DEPENDANCE 

The dependance of the local age parameter s(r) on 
lateral distance exhibits the typical dependance at the 
level of Akeno experiment with a minimum near 30-40 
m distance and larger values near at small radial dis- 
tances or at large distances from the axis. An example 



FIG. 3: Lateral electron distributions at Kascade level. 
Experimental(full line), NKG option(dashed), NKG- 
Lagutin modified function dotted 

is given in the fig. 4 for N e = 10 6 . 

The comprehensive dependance of the lateral age 
parameter s on N e is shown on fig. 5 

The lateral age parameter is obtained as the value 
of the La. p. s(r) averaged inside the shower disc up 
to the shower radius Ru m (distance where the average 
lateral distribution reaches a threshold density satis- 
fying the trigger conditions and remains detectable by 
the scintillator detector used i.e. for instance lel/m 2 
for a counter of lm 2 area). One of the most simple 
interpretations of the minimum near the knee on fig. 5 
could be a mixed component progressively enriched in 
nuclei after the knee as suggested half a century ago 
by the galactic leakage, according to the different Lar- 
mor radius of nucleons and heavier nuclei. However, 
for small energies the shower radius Ru m is lower than 
40m whereas it exceeds 100m at large energies above 
the knee; the calculation of the average s± up to those 
different limits may be the reason for decreasing s be- 
low the knee and an increase in age beyond the knee. 



IV. CONCLUSION 

The experimental profiles cannot be expressed in 
terms of NKG function with a single age parame- 
ter. Consequently, though the lateral age parameter 
contains information of both hadronic cascading and 
primary composition, its determination can be biased 
while fitting the experimental distributions. The char- 
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Local age parameter in Akeno 



acteristic dependance of local age parameter with ra- 
dial distance allows a more accurate determination of 
s± as well as the conversion to sl 
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A detailed study of dependence of s± on shower size, 
included in a multiparameter analysis (muon-electron 
abundance, absorption length) and comparison with 
the Akeno and the Kascade data (in progress) should 
help to understand the knee region in terms of primary 
composition. 
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FIG. 4: Local age parameter versus distance at Akeno. 
Experimental points squares (thick scintillators), stars 
(thin scintillators), dotted line proton primaries, dashed 
line iron primary, horizontal lines longitudinal age param- 
eter respectively for p and iron 
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FIG. 5: Lateral age parameter versus size N e at Akeno 
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